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We directly combine ideas of the quasiclassical approxi- 
mation with random matrix theory and apply them to the 
study of the spectrum, in particular to the two-level corre- 
lator. Bogomolny's transfer operator T, quasiclassically an 
N x N unitary matrix, is considered to be a random matrix. 
Rather than rejecting all knowledge of the system, except for 
its symmetry, [as with Dyson's circular unitary ensemble], 
we choose an ensemble which incorporates the knowledge of 
the shortest periodic orbits, the prime quasiclassical informa- 
tion bearing on the spectrum. The results largely agree with 
expectations but contain novel features differing from other 
recent theories. 



The study of wave systems of low symmetry is known, 
somewhat misleadingly, as "quantum chaos" . Aside from 
experiment and numerics, the most important tools in 
this field are the quasiclassical approximation [QCA] [jjj , 
and random matrix theory [RMT] 

These two approaches are complementary. Standard 
RMT treated problems where knowledge of the Hamilto- 
nian is almost nil. It identifies universal features, notably 
of the spectral statistics, which statistically persist for all 
quantum systems with time reversal invariance, say. 

Standard RMT fails for the most fundamental spectral 
statistic, however, namely the mean or smoothed [Weyl] 
density of states, d(E), which is not, after all, so broadly 
universal. However, it is argued that local quantities can 
be calculated, i.e. quantities involving energies in the 
range E ± W, where W d(E) » 1, and W is small 
enough so that d{E) can be taken as constant. RMT 
treats the Hamiltonian, TL, as a random N x N matrix, 
locally normalizes to d, and lets N — > oo. 

In the QCA, a system is chosen, e.g. a particular 
stadium billiard, and the smallness of the ratio [wave- 
length] /[characteristic classical length] is exploited. A 
particular spectrum, not just statistics, is in principle 
obtained. The basic idea, of Gutzwiller j3j but much ad- 
vanced by others Q] , || , |(| , , makes two subtly related 
but ostensibly different QCA approximations. 

The first is the Weyl approximation to d whose 
leading order is the phase space integral d(E) = 
(2n%)f I d*fpdfq5(E — H(p,q)) where the dimension / 
is here taken to be 2. The density of states is divided 
into d(E) = J26(E - E a ) = d(E) + d osc (E). A 
WKB-like approximation is made to give d osc , describ- 
ing the oscillations about the Weyl term, d osc (E) — 
A p cos S p (E)/h. The sum is over periodic orbits of ac- 



tion S p [which includes the Maslov index] and A p depends 
on the classical stability of the orbit. In this primitive 
form, the result is the Gutzwiller trace formula jj| . 

Many advances have been made, e.g. d osc can be writ- 
ten as a ratio of absolutely convergent series, rather than 
as the divergent series above The major result is that 
by such means only orbits whose period T p = dS s /dE is 
less than half the "Heisenberg time", t% — 2nhd oc 1/Ti 
need be considered: All longer orbits can be effectively 
expressed in terms of these. Orbits longer than com- 
bine to give nothing, those shorter display a sort of re- 
flection symmetry [resurgence] about ^Tn- 

No direct way to solve a fundamental problem is 
known, however. Namely, the number of periodic orbits 
proliferates exponentially with T p and for small ti no one 
can actually find such long orbits, let alone add them up. 
Therefore recourse is made to a statistical treatment of 
the long orbits, via the Hannay-Ozorio de Almeida sum 
rule ||, valid for hard chaotic systems. This requires a 
form of the diagonal approximation [DA]. [The product 
doscdosc is expanded in a double Gutzwiller series, and 
only the diagonal terms are kept. After the DA, quan- 
tum effects enter only trivially through the energy scale 
d] The DA by itself is not sufficient, and it fails dra- 
matically [at long time] since long orbits need to show 
important correlations, if they are to describe a discrete 
spectrum, and give results in agreement with RMT and 
experiment for neighboring level properties pQ|. [ A re- 
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cent effort of Bogomolny and Keating 
discreteness of the spectrum in an ad 
way. The supersymmetry [SS] technique has begun to 
find good results in this connection also ]l^| .] In short, 
the QCA has to invoke statistics to deal with the long 
orbits. On the other hand, there are many influences 
of the short periodic orbits on the spectrum which are 
observable and have been of great interest. 

A crucial advance in semiclassics was made by Bogo- 
molny H with his T operator. The T operator is an 
approximation to and/or generalization of the boundary 
integral method used to study billiard problems. It is a 
QCA to an exact kernel K(q,q',E) [where q,q' run over 
the boundary of the billiard] such that the spectrum is 
given by the zeroes of det(l — K(E)). The equivalence 
of the 'resurgence' mentioned above was shown: Namely, 
in QCA T is unitary and of rank N , i.e. equivalent to 
an N x N unitary matrix T mn (E), Further, a general 
Poincare surface of section can be used instead of the 
boundary, extending the technique beyond billiards. 

Our proposal is to treat the matrix T mn as a ran- 
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dom unitary matrix, generalizing Dyson's famous work 
H . This has several advantages over the usual choice of 
the Hamiltonian as random matrix. 

• The mean distribution of eigenvalues e l6n , d(6), is 
physically meaningful, and related to d{E) . [We dis- 
tinguish the distributions in angle and energy by 
the argument. We also denote, by d without an ar- 
gument, the Weyl density of states. Dyson's COE 
and CUE matrix ensembles have d(9) = N/2%.] 
The study of deviations from d constant is the pur- 
pose of d osc . Random Hamiltonian matrices are 
unable to deal with such questions. 

• The rank N is meaningful. It is [usually] given by 
N ~ 2L/X(E) where L is the length of the surface 
of section, and A is the wavelength at energy E. 

• The matrix size is N oc h~~ 1 rather than h~~ 2 , an im- 
portant efficiency of the boundary integral method. 

• Specific quasiclassical knowledge may be incorpo- 
rated into the definition of the ensemble of matri- 
ces in a natural and convenient way. In particu- 
lar, one may study modified Dysonian ensembles, 
which have prescribed mean values of the traces 
oy = TrT r . The set of traces oy ; r = 1..N uniquely 
determine det(l — T) and thus the spectrum. Low 
order traces are usually well evaluated quasiclas- 
sically by the method of stationary phase. They 
are expressed in terms of short periodic orbits, and 
we take this approximation for oy. The high order 
traces cannot usually be so evaluated as the number 
of orbits is too large. 

Some difficulties with the approach are 

• Direct calculation gives the eigenphases 9 n rather 
than the energies E a for which 6 n (E a ) — mod 2tt. 
However, there is evidence that the 9 n are locally 
linear functions of the energy. For N large, this 
is sufficient, as Dyson showed. [Bogomolny, in un- 
published lectures, has made this connection more 
precise.] For billiards one can probably do better, 
since scaling indicates that n = a(k — k n ) where 
k = s/E. [We put h = 1.] 

• We have assumed that N is an integer which de- 
pends on energy. Therefore, rather unphysical 
jumps in TV are needed as the energy increases. Per- 
haps the simplest way around this is to express the 
results as a function of N and 'analytically con- 
tinue' to continuous N. 

Let Z{E) = det(l - T) = - e 4e "). Then 
d osc = #Im^ \nZ(E) = -J-J2°'n( 1 + Imcot ¥n) 

n 

= z2i^'n+ d 'J(UE))} (1) 



Assuming 9 n (E) = (E - E n )9' n with 9' n _ « 9' ap- 
proximately independent of n, gives d = dweyi(E) = 
SFErX = If 61 '- For a billiard, Q< = a/2k. Using 
the estimate for N it is found that d = La /An 2 . Thus 
a = Att/L, where A is the area of the billiard. 

We now add knowledge of the short period orbits 
and only this knowledge to the problem, in the spirit 
of the information theoretical approach of Balian [13 
We cite other recent work using this technique [14 
Let dfi(T) be the measure for the class of random ma- 
trices under consideration. We want the probability 
distribution P(T)d(i(T). Then the 'entropy' S[P] = 
- J dfi(T)P(T)ln(P(T)) is maximized, subject to con- 
ditions / P{T) TrT r dfi(T) = oy; r = 1, 2, .. 

For simplicity we choose the measure without time re- 
versal symmetry corresponding to Dyson's CUE. T can 
be parametrized by its N eigenphases 9 n together with 
N(N — 1) other real parameters b s , which we integrate 
out. Then d(j,(T) = U„< m ~ e w ^\[ n dB n . Intro- 
ducing Lagrange multipliers for the constraints, we ob- 
tain 



p(9i,9 2l ..9 N )l[d8 n = J] |e*' -e*-| 8 JJ«j(0„)d0 n 



where the positive weight factor is 

w(9) — exp X r cos r(9 — ?9 r ) 



(2) 



(3) 



The A's and i?'s are the Lagrange multiplier parameters 
[two for each complex oy] needed to specify the traces 
and A r vanishes for an unconstrained trace. We study 
d(9), R(9i, 6*2) obtained by integrating out all but one or 
two 0's in Eq.(§). 

The method of orthogonal polynomials is used [0, 
|^5| . Polynomials 0„(z)in z = e 10 are found, satis- 
fying ^ / dOw (9) 4> n (z) 4> m {z) = 8 nm where the bar 
indicates complex conjugation. This and subsequent 
integrals are over [— tt, tt]. Eq.(||) can be expressed 
p(9 ll ..9 N ) = detM+M = detK where V2nM nm = 
$„(M - v^Mi(^) and K H = K(8 t ,9 J ) = 
If Yln=o ®n(9i)® n (9j) . The reason for writing p in this 
form is that an integral over one variable of an R x R 
determinant of the form det K is proportional to the 
{R — 1) X [R — 1) determinant over the same argument. 
We thus find that 



and 



d{9)=K(9,9) 



R(6 1 ,6 2 ) = d(6 1 )d(6 2 )-\K(6 1 ,6 2 )\< 



(4) 



(5) 



The sum in the definition of K can be carried out [Jf5[ 
The formula [Christoffcl-Darboux] is 
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4>n(zi)4>n(z2) 



z{ 4>n(z 1 )z 2 </>n{z 2 )-4>n{zi)<Pn(z 2 ) 



1 - z x z 2 



(6) 



[4>{z) = 4>{z)] We thus express everything in terms of the 
single polynomial </>jv [and w). 

Given w(9), satisfying certain conditions, a formula ex- 
ists for 0jv for N — > co. The weight function of Eq.(|j|) 
is particularly congenial. If the exponent is a convergent 
series, the conditions are well satisfied. 

The weight can be uniquely expressed by w(9) = 
\D(z)\ . Ambiguity in D is removed by insisting that 
D{z) for \z\ < 1, and D(0) = 1. For Eq.(|), 
D(z) = exp [| \ n e~ m§n z n ~\ ■ Asymptotically, for large 
N, the theorem 11151 is that 



[D(z 
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Assuming asymptotia has been reached, we use Eq.(0) in 
Eq.(0) and Eq.ffl) to find, with L'Hopital's rule, that 



d(9) 



N 
2^ 



Re 



D'{z) 



D(z) 



(8) 



The integral over the nontrivial term vanishes. Thus 



d{9) 



N 
2^ 



-y 



nX n cosn(9 — i3 r 



(9) 



The trace conditions are ay = Jd9e lr9 d(9) = 
^rX r e lrlSr , determining the Lagrange parameters. The 
mean density of states in energy is d(E) = 9'd(9 = 0) = 
d+ — \&r\ cos(rd r ). Gutzwiller's density of states, ex- 
pressed in terms of the traces of T is 



d Gutz (E) = d+ ilm^^ia r (^) 



(10) 



where for pedagogical purposes we assume that e r is 
rapidly varying. If the derivative d' r ss 9' , the leading 
terms of dcutz coincide with the RMT d(E). 

We thus have the interpretation that, at the level of the 
mean density, the random T operator scheme is equiva- 
lent to moving the long wavelength part of d osc into the 
mean density of states. The formulas above give a way 
of including this knowledge into a prediction of the cor- 
relation functions. 

We next evaluate the two point correlations, using 
Eq.(g). Define ip{d)[= \ £ A„ sinn(6> - d n )\ by e^W = 
D(z)/ \D(z)\ . [For Dyson, tp = 0.] Then 

- f/)(02) = y K sin f (0i - 02) cos § (0! +9 2 - 20„) 



\K(9 1 ,9 2 )\ 2 = 



[%(0 1 -0 2 )+1>{0 1 )-il>(02))] 



2w sin i ( 



-62) 



(11) 



agreeing with Dyson if if) = 0. We wish to calculate 
C(x) = (C e (x)) w = ((^(e+ £)d(e- 

= ((?E^+s-w-S-^ 



a j 6 




(12) 



The outer average is over the energy range E ± W. The 
inner average is over the random matrix ensemble. 
We assume as above that 9 a (E) = (E — E a )9' . Then 



R(E U E 2 ) = {^) 2 R((E-E l )9',(E-E 2 )9') 
Using this in Eq.(n2[), with e = E we obtain 



(13) 



C ^=d s{ d )+ {1) d {--d) d \-d 



e'Nx 



2d 



J2K sin^ff cosm? 



2nd , 



2d 



(14) 



Since 9' /d — 2tt/N is small, we find (suppressing S(x)) 

c e (x) = (i+5> n+ ) (i+y^ 

sin ttx (l -I- J2 !J Ji a - cos n ^n 



Tr.r 



(15) 



where c n ±(x) = ^cosn(^f ± n ) 

We now perform the energy average. Assume that 
only the $ n 's vary significantly, and that they can be 
averaged independently. This gives the DA for the 
squared terms c„ + c„_ in the first line of Eq.jT^), equiv- 
alent to that found by Berry jlf}]. We call this term 

Cb(x) = |E(^) 2 cos2fix. Let p(y) = p(-y) = 
{$(y — S cosn $n)) w be the distribution of the cor- 
rection term in the second line of Eq.(p^|). Its width is of 
order A -1 . The form factor |l(| is 

k(t) = l + J dxe 2 ^C(x) = ^2{^) 2 S(\t\-tn) 



t\-i 



dy / dy'p(y' 



(16) 



where t n (= n/N) is the period of the orbits associated 
with the n'th trace. The first term is a simulation of 
Berry's DA result for the short orbits. The second term 
rounds off the sharp corner found in standard RMT, 
[where p = S(y).] The magnitude and width of the round- 
off are both ~ A -1 . Since p > 0, we do not find the 
suggested Q oscillations in k(t) near t = 1. 

This expression mostly agrees with what is believed to 
be the correct form of k(t). For large \t\, the discreteness 
of the spectrum dominates. For \t\ ~ t c ~ A -1 , the DA 
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applied to Gutzwiller's expression gives spikes when \t\ is 
at a period of a short periodic orbit. We should improve 
our elementary estimate to replace Cb as given above by 
Berry's result, and we should also make the correspond- 
ing change in ip(9). For t c < t « 1, the DA applies 
and the sum rule |^] gives a straight line of unit slope. 
There is then structure near \t\ — 1 which is a conse- 
quence of the short period structure. The area removed 
under k(t) near t = 1 just compensates the area gained 
from the spikes at t ~ t Cl a consequence of R(9, 9) = 0. 
It's unclear whether resurgence plays a role in produc- 
ing structure near \t\ = 1, as has been suggested [ jl2| . 
Our formulation contains all the ingredients needed for 
resurgence. RMT fluctuations in det T may wash out this 
effect. It seems rather that the energy average produces 
its biggest effect at |t| = 1 because of the kink there. 

The replacement sin mix /N — > nnx/N in obtaining 
Eq.(|i~5|) from Eq.(14) is not necessary. Averaging as be- 
fore, we find 



/ l n \ ^ . irnx 

^cos \ 2kx + 2 A n sin — — cos i 



n( 



cos 2ttx I I ( 1 — Ad sin 



2 ^) =x 2 V(x)cos2irx (17) 



and we have assumed A„ small. Then 
1 / 1 <9 2 lnX>(x) 



2tt 2 V 2 dx 2 



+ V{x) cos 2irx , (18) 



the relation found by SS |2|. With Eq.(|17|) small struc- 
tures near \t\ = 1 decorate the result of Eq.(|I^) but the 
result is still monotonic. 

It is not clear under what conditions Eq.(|l^) holds. 
The SS result is derived for large x, while ours is for 
small A„, in other words, for rather unstable orbits. 
The V{x) corresponding to a distribution p, is T>(x) — 
x~ 2 J dy p(y) cos 2irxy, but the term with h\T> does not 
reproduce Cb- The T> of jl2| is derived from Cb as 
a correction to the sum rule: it corresponds to a non- 
positive 'distribution' p(y) oc cos 2^722/ exp (—27r7i |y|) 
where 71 + 172 is the lowest nonvanishing eigenvalue of 
the Perron- Frobenius operator. It of course does not give 
the spiky structure of k(t) at t ~ t c . 

To summarize, an application of random matrix theory 
to Bogomolny's T operator, one which retains knowledge 
of the short periodic orbits, leads to a relatively simple 
theory which gives reasonable results for the density of 
states and for the two point correlations. It does not 
agree in every detail with previously obtained results, 
however. This may be due to approximations we have 
made as well as the fact we have not specifically consid- 
ered orbits of length of order 7-f 1 > t c . 

In this note, we have not worked out what happens 
for finite N, although the results should be meaningful. 
One can find the results for higher correlations directly. 
We don't expect them, or the corrections to the nearest 
neighbor level spacing, to be very illuminating, however. 
We hope to extend the results to other symmetry classes, 



and to parametric correlations. It should also be possible 
to compare statistics of eigenphases directly with this 
theory, in the case of quantum maps, [and indeed the T 
operator is the quantum Poincare map.] 
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